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Gamow- Teller(GT) transitions in the (5 decays of ^Li and ^^Li are investigated with theoretical 
calculations of antisymmetrized molecular dynamics. The calculated B{GT) values are small for 
transitions to low-lying states in Be isotopes, while relatively large B{GT) values are found for 
excited states at excitation energy Ex > 10 MeV. Sum of the B{GT) values are discussed for each 
spin parity of final states. The calculated results seem to be inconsistent with the experimental 
report of the strongest GT-transition from ^Li to the 5/2~ state at 11.8 MeV of ^Be. 
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I. INTRODUCTION 
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Experimental studies of (3 decays from unstable nuclei near the drip lines are developing nowaday due to recent 
progress of experimental techniques. Gamow- Teller(GT) transition strength B{GT) values were extracted not only 
for transitions to low-lying states but also those to highly excited states. Moreover, determination of the B{GT) 
values has been performed also in high resolution experiments of charge exchange reactions as well as the /3 decays. 

In these years, measurements of the (i decays of ^Li and ^^Li were carried out in several experiments, for example, 
at the ISOLDE facility in CERN. They provide information of new states in ^Be and ^^Be such as the excitation 
energies, spins and B(GT) values. In the (3 decay of ^Li, the GT transitions to low-lying ^Be states in < 10 MeV 
are weak while the strong transitions to the 11.81 MeV state were re ported[l|, i]. The extracted B{GT) value for 
the ^Be state at = 11.81 MeV is surprizingly large as B{GT) — 8.9(1.9) which takes 65 % of the Ikeda sum rule, 
'^{qa/ gvY' ~ 13.4. In a new measurement of the /3 decayQ, its spin and parity were assigned to be 5/2~. Compared 
^ with the mirror transitions from ^C to the ^B(5/2~) state at 12.19 MeV;4j, this suggests the abnormally large mirror 
Ch . asymmetry that the B{GT) value for the ^Be(5/2~) at 11.81 MeV is larger by factor 4 than that for the mirror 
transition as argued in Ref. [sj. However, the shell- model calculations in ap-shell basis do not indicate such the large 
B{GT) values in the ^Li decays nor the large mirror asymmetry |5i|. 
^ ■ The experimental measurements of the j3 decay were performed also for the drip-line nucleus, ^^Li[l, 0, Many 
l/^ \ excited states of ^^Be were observed in low energy region. The GT transitions to ^^Be states in E^ < 10 MeV 
are not strong and the extracted B{GT) values are less than B{GT) = 0.5 for the low-lying statesQ. In a recent 
' measurement, a new state at 18 MeV in ^^Be was observed and the value B{GT) = 1.9 for this state was extracted 
. . from the branching ratiosQ. 
0^ ■ As mentioned above, in the (3 decays of neutron-rich Li isotopes, the GT transitions to low-lying states of ^Be and 
^^Be are relatively weak compared with rather strong GT transitions to highly excited states. It is understood easily 
because of the cluster features of Be isotopes. In theoretical studies of Be isotopes, it was suggested that 2a-\-n and 
2a-f 3n cluster structures develop in most of the low-lying states of ^Be and ^^Be, respectively [lol [Til [T2|. [l3l [l3. [Tsl. [l^ . 
The GT transitions to such the cluster states should be suppressed because GT transition to an ideal a cluster written 
by the (Os)^ configuration is forbidden exactly due to Pauli principle. In other words, the high-lying states in E^ > 10 
MeV with the strong GT transitions are expected to have significant components of cluster breaking. 
' In this paper, we investigate the GT transition strengths of the /3 decays, ^Li^^Be and ^^Li— >^^Be, with a theoretical 
method of antisymmetrized molecular dvnamicsfAMD) [TtI. [isl [l9| . Distribution of the B{GT) values is analyzed for 
each spin of the final states. Particular attention is paid to the strong GT transitions to the ^Be 5/2— state at 11.81 
MeV suggested in the experimental report. 

This paper is organized as follows. In the next section, we describe the formulation of the present calculations. We 
explain the adopted effective interactions in lllli and show the calculated results as well as the experimental data in 
IIVI Finally, in |Vl we give a summary and an outlook. 



o 



II. FORMULATION 



We use a method of AMD which has been proved to be one of the useful approaches for structure study of light 
unstable nuclei. In the present work, we first perform the variation after projection(VAP) with respect to the total- 
angular-momentum and parity projection(spin-parity projection) in the AMD and extend the method to calculate 
GT strength functions to final states in a wide e nergy re gion. The detailed formulation of the basic AMD method for 
nuclear structure study is described in Refs. p^[l8lTl9t. 
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In the AMD method, a wave function of a j4-nucleon system is written by a Slater determinant of single-particle 

Gaussian wave packets, 

$amd(Z) = -^A{ipi,ip2, ^Pa}, (1) 

where the ith single-particle wave function is written by a product of spatial((/)), intrinsic spin(x) and isospin(r) wave 
functions, 

= (h^iXin, (2) 

<?i'x.(rj) cx exp{-i/(rj- - — (3) 

Xi = {^+^i)x^ + {\-ii)Xl- (4) 

(j) and X arc represented by complex variational parameters, Xii, X2i, Xa^, and The isospin function is fixed 
to be up(proton) or down(neutron). We take a fixed width parameter u which is optimized for each nucleus. That 
\s V = 0.20 fm~^ for ^Be and ®Li, and v = 0.18 fm~^ for ^^Be and ^^Li in the present calculations. Accordingly, an 
AMD wave function is expressed by a set of variational parameters, Z = {Xi, X2, • • • , X^i, ^1, C2, • • • , £,a}- 

For the lowest J"^ state, we vary the parameters Xj and = I ^ A) to minimize the energy expectation value of 
the Hamiltonian, with respect to the spin-parity projected AMD wave function; $ = Pj^'5^$AMD(Z). 

Here, P//^ is the spin-parity projection operator. After the energy variation, the optimized parameter set Zf^ of the 
minimum-energy solution for the lowest J"^ state is obtained. The solution Z;^'^ for the nth J'^ state is determined by 
varying Z so as to minimize the energy of the wave function orthogonalized to the lower states; 

n-l 

1$) = |P^'^$amd(Z)) - J2 l*i^^^")(*i^^^n^i?^K^AMD(Z)), (5) 

fc=l 

where ^^^^^ normalized wave function determined for the lower states. This is the standard procedure of the 

VAP calculation in the AMD method. 

In the VAP calculations, the wave functions for the states are obtained one by one from the lower energy states 
and it is not easy to calculate all final states of the GT transitions which fragment generally into many high-lying 
states. In order to exhaust the GT strengths from an initial state, we extend the basis wave functions by operating 
one-body spin-isospin operators to the VAP wave functions for the parent nucleus. Let us consider the f3~ decay of 
^Li. We perform the VAP calculation of the lowest J'^ = 3/2~ state of ^Li and obtain the wave function ^amd(^i^^~) 
for the ^Li ground state. Here we rewrite 7?^'^~ for ^Li with Zi„it. Final states of the GT transitions from the ^Li 
ground state are = l/2~,3/2~ and 5/2^ states in the daughter nucleus '^Bc. As for the final states in ''Be, we 
prepare the basis wave functions in two steps as follows. We first do the VAP calculations for the lowest two 
states of ^Be, and obtain the wave functions $amd(Z;!'') for = 1/2J", 1/2^, 3/2j", 3/2^, 5/2j" and 5/2^ which 
approximately describe the corresponding low-lying states. Next, we create other basis wave functions by operating 
the one-body spin-isospin operators to the ^AMoC^init) as, 

*fc,a(Zinit) = ■^^-4{(pi,</92,---,0-aT"<y5fc,---,</9A}, (6) 

where </5fe (fc = 1, • • • ,6) is a neutron wave function in ^Li. ctqT" is the one-body GT transition operator where 
Ga.[pi = x,y and z) is the spin operator and t~ is the charge changing operator. As a result, we get 6 x 3 = 18 number 
of basis wave functions, ^^'^(Zinit), which exhaust the GT transition strengths from the ^A^MoCZinit)- 

Finally we determine the wave functions of the final ®Be states by the spin-parity projection and the superposition 

of all the 6-1-18 basis wave functions, 

^Jn = E <^^J'%rnPMK'^'lMu{'Z'in + E c£^^MK$fe:«(Zi„H). (7) 
J' Km Kka 

The coefficients c^ji^]^ „ and c^^^^ are determined by diagonalizing Hamiltonian and norm matrices. 
The value B{GT) of the GT transition strength is written by the square of the reduced matrix element 

B{GT) = {gA/gvmof^-^)\'' (8) 
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of the GT transition operator, 

(yiGT±) (9) 
i 

where gA and gv are the axial- vector and vector coupling constants and taken to be {gA/gv)'^ = 1-51 in the present 
calculations. It is worth mentioning again that the GT transition strengths from the ^Li ground state written by the 
wave function Pmk^^^~ ^^AMnC^i^ii) exhausted by the final states, <l>j„ (n = 1, • • • , 24), of ^Be. 

III. EFFECTIVE NUCLEAR FORCES 

The effective nuclear interaction adopted in the present work consists of the central force, the spin-orbit force and 
the Coulomb force. For the central force, the MVl force case (3) [2l[ containing finite-range two-body and zero-range 
three-body terms is used. As for the spin-orbit force, the same two-range Gaussian form as that in the G3RS force 
[121 is adopted. 

The used values of the interaction parameters in the MVl force are b = h — and m — 0.62 of Bartlett, Heisenberg 
and Majorana parameters, and the strengths of the spin-orbit force are taken to be uj — —uu = 3000 MeV. These 
parameters are the same as those used in the AMD4-VAP calculations of ^^C and ^°Be in Refs. ^^[2^. Hereafter, we 
call this parameter set (A). In order to see the interaction dependence of the GT transition strengths, we also use two 
sets of modified parameters, (B) and (C). The set (B) b = h = 0, m = 0.62 and ui = —uu — 2000 MeV has the same 
central force but the weaker spin-orbit force than (A), and the set (C) b = h = 0.15, rn = 0.62 and uj = —ujj = 3000 
MeV contains the Bartlett and Heisenberg terms in the central force but the same strengths of the spin-orbit force as 
(A). 

IV. RESULTS 
A. GT transition from ''Li to ''Be 

In the present results of the VAP calculations for low-lying states of ^Be, it is found that 2a + n cluster structure 
developes in the ground state(3/2^) and excited states, 5/2^, l/2j", 3/2^ and 5/2^. The calculated energy spectra 
of these low- lying states are consistent with those calculated with the 2a + n cluster models [l3, [lil, [3 and the 
experimental energy levels (table |T|. 

By adding the basis wave functions <I>^'^(Zinit) constructed from the spin-isospin operated ^Li wave functions to 
the VAP wave functions as explained in|TTl the final wave functions of ^Be are obtained and the B{GT) values for 
the transitions ^Li^^Be are calculated. The calculated B{GT) values to the final ^Be states up to ~ 13 MeV are 
shown in tableU The B{GT) values to the ground state(3/2j~), the 5/2j~, l/2j~, 3/2^ and 5/2^ states are small. This 
is because those states of ^Be have the 2a -\- n cluster structure and they have small overlap with the initial state of 
^Li having no developed cluster structure. Fig. [T] shows the density distribution in the intrinsic wave functions of the 
^Be ground state(<I'^^j3(Z^^^~)) and that of the ^Li ground state(<I'2MD(^iiiit))- It shows that the intrinsic structure 
changes drastically from the initial state to the final state. Moreover, the GT transitions to such the cluster states 
with the 2a core are suppressed because of Pauli principle as mentioned before. In other words, the finite B{GT) 
values account for the slight dissociation of the ideal a clusters in ^Be. 

In contrast to the small B{GT) values to the low-lying states, the excited states at ~ 12 ^ 13 MeV show 
significant GT strengths which come from the non-a-cluster states constructed from the spin-isospin operated ^Li 
wave functions given in Eq. [6l The strongest GT transition is found for the 1/2^ state in the present results. In the 
experimental measurements of the /S decay of ^Li, rather large values of B{GT) were reported for the states around 
E^ = 12 MeVP, [2, [1 as B{GT) = 1.44 for the E^ = 11.28 MeV state and B{GT) = 8.9 for the E^ = 11.81 MeV 
state. The former state is suggested to be a 3/2 ~ state [lij and it might correspond to the theoretical 3/2;^ state in 
the present result. The calculated value of B{GT) = 0.94 agrees to the experimental value. The spin and parity of 
the latter E^ — 11.81 MeV state was assigned as J'^ = 5/2~ in the recent analysis of the (3~ decay of ^Li[3|. Although 
this state seems to correspond energetically to the calculated 5/2;^ state, the theoretical B{GT) value is much smaller 
than the experimental value of B{GT) — 8.9 which exhausts 65 % of the Ikeda sum rule. As shown later, there is no 
5/2~ state having such the strong GT strength compatible to B{GT) = 8.9 in the present calculation. Also in the 
shell-model calculations, it is difficult to reproduce the extraordinary large B(GT) value of the experimental dataQ. 

Figure [5] shows the B{GT) distribution to excited states of ^Be up to E^ = 30 MeV calculated by using the 
interaction set (A). In the B{GT) values to l/2~ states, the largest peak is found at E^ ~ 12 MeV, and some 
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strengths distribute around 20 MeV. The transition strengths to 3/2~ states distribute mainly in = 10 — 20 MeV 
region. The B{GT) values to 5/2~ states are relatively small compared with those to l/2~ and 3/2~ states. In 
order to show the interaction dependence of the GT strengtions, we show the B{GT) distribution calculated with 
three interaction sets (A), (B) and (C) in Fig. [31 The B{GT) distribution is qualitatively similar between these three 
interaction sets, though the broadness of the distribution changes slightly. In all the results, the B{GT) values are 
very small for the transitions to low-energy Ex < 10 MeV states, while the strengths distribute mostly in Ex = 10 — 20 
MeV. Compared with the shell model calculations in Ref. 25], which show the B{GT) distribution concentrating at 
Ex = 10 — 12 MeV, the present calculations show significant fraction of the B(GT) values to highly excited states in 
Ex>lh MeV. 



decays of C and Li, and theoretical values of the Li decay. The 



TABLE I: Experimental data of the GT strengths in the 

experimental data of the B(GT) values are taken from [2l and ''[4|. "The spin and parity 3/2~ of the ^Li(11.28 MeV) is the 
assignment of Ref. [241 but it is not confirmed yet and other spin-parity candidates, 7/2~ and 7/2 , are sug gestedHi]. "=5/2- 
was assigned in Ref. [3|- 
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FIG. 1: The distribution of matter, proton and neutron density in the intrinsic wave functions of (a)<I>XMD(Zi'''^ 
ground state and (b)$AMD(Zinit) for the ^Li ground state calculated with the interaction (A). 



for the Be 



B. GT transition from "Li to "Be 



The B{GT) distribution in the ^^Li^^Be decay calculated with the interaction (A) is shown in Fig.[31 The B[GT) 
values to 1/2" states concentrate around Ex = — 18 MeV, and those to 3/2" states distribute widely in various 
excited states. The results of the B{GT) distribution calculated with three sets (A), (B) and (C) of the interaction 
parameters are displayed in Figs. [Hl^a), [Sj^b) and[n]Jc) as well as the observed B{GT) values in Fig.[ni[d). The results 
are qualitatively similar between these three interaction sets. The B{GT) values are small in Ex < 10 MeV, while 
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FIG. 2: The B{GT) distribution to (a)l/2" states, (b)3/2~ states and (c)5/2" states of ^Be in the P decay of ^Li calculated 
with the interaction (A). 



they distribute widely in = 10 — 25 MeV region and shows the broad peak structure with the center around 
Ex = 15 — 20 MeV. The small B{GT) values to low-lying states are understood by the 2a core structure in ^^Be, 
while the highly excited ^^Be states with significant GT strengths are those with non-cluster or less-cluster structures. 
This is a similar situation to that in ^Be. 

Although ^^Li is known to be a neutron-rich nuclei with the neutron-halo structure, halo efTects are not taken into 
account in the present calculations because single-particle wave functions in the present model are restricted to be a 
Gaussian form and are not suitable to explain the long tail of the halo-neutron wave functions. The halo effects in 
the B{GT) values for the "Li-*"Be decay were discussed in the shell model study of Ref . [l^l , which showed that 
the halo structure gives small quenching of the B{GT) values. 



C. Sum rule 



The Ikeda sum rule for the GT transition strengths is given as, 

SiGT^) _ SiGT+) ^ 4(0|[O(^^+',O(^^-^]|0) = i^i{N~Z). (10) 
9v 9v 

Here S^^^^'^ is the sum of B{GT) values for the P'^ decay and |0) is the initial state. Since 

5'(GT+) fQj, ^Yi^ p+ decays 

of neutron-rich nuclei is very small, the sum rule is often approximated as S'''^"^^-' ~ ^3{N — Z) by neglecting the 

5'(GT+) tgj.^^^ 

The values of the sum rule 3(A^ — Z) x (gA/gv)^ are 13.59 and 22.65 for the ^Li and ^"'^Li f3~ decays, respectively. 
The theoretical values of sum of the B{GT) values to l/2~, 3/2~ and 5/2~ states in Be are listed in table [TTl It is 
important that the sum values are determined only by the initial state of the Li ground state but do not depend on 
final states of Be. Compared with the results (A), (B) and (C), it is found that interaction dependence of the sum 
values is small. It means that the wave functions of the ^Li and ^^Li ground states are not sensitive to the adopted 
interactions in the present calculations. The sum of the B{GT) values to l/2~ states is the largest. On the other 
hand, that to 5/2" states is the smallest and exhausts only one forth of the calculated total sum S*'*^-^") value in both 
the ^Li and ^^Li decays. This result contradicts to the experimental report of the value B{GT) — 8.9 to the 5/2" 
state at E^ = 11.81 MeV in Ref. Theoretically, the sum of the B{GT) values for each spin is determined by the 
structure of the ^Li ground state with the spin and parity 3/2", and it is unnatural that the highest-spin 5/2" state 
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FIG. 3: The B{GT) distribution for the decay ^Li^^'Be. The B{GT) values to 1/2", 3/2" and 5/2" states of ^Be are piled 
up. The upper panels (a),(b) and (c) show the results calculated with the interaction sets, (A), (B) and (C), respectively. The 
lowest panel (d) shows the experimental data which are taken from the same references as those in table |l] 



exhausts such the large fraction of the Ikcda sum rule. In other words, the experimental value B{GT) = 8.9 for the 
5/2~ state seems to be too large to be described by theoretical calculations if the ^Li ground state has an ordinary 
structure. 

TABLE II: The calculated results for sum of the B{GT) values of ^Li and ^^Li /? decays to J~ states in ^Be and ^^Be. The 
values of Ikeda sum rule 3(A'^ ~ Z) x [gA/gv)'^ are 13.59 and 22.65 for ^Li(/3^) and ^^Li(/3~), respectively. 
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FIG. 4: The B (GT) distribution to (a) 1/2 states, (b) 3/2 states and (c)5/2 states in the ^^Li^^^Be decay calculated with 
the interaction set (A). 



V. SUMMARY AND OUTLOOK 



The Gamow-Teller(GT) transitions in the /3~ decays, ^Li^^Be and ^^Li-^^^Be, were investigated with a method 
of antisymmetrized molecular dynamics. The calculated B{GT) values are small for transitions to low-lying states 
of Be isotopes because of the 2Q;-core structures in the final states. Significant strengths are found in the B{GT) 
distribution to non-cluster states of ^Be in > 10 MeV region. Sum of the B{GT) values for each spin parity of the 
inal states was also studied. 

Particular attention was paid to the strong /^-transition from ^Li to the ^Be(5/2^) state at 11.8 MeV which was 
reported in the experimental workQ. The present results are inconsistent with the strong GT transition to the 5/2" 
state which shows the large fraction of the Ikeda sum rule value. Also in terms of the sum rule, the experimental value 
B{GT) = 8.9 for the 5/2" state seems to be too large to be described by theoretical calculations if the ^Li ground 
state has an ordinary structure. 

This work is the first achievement in which the AMD method was applied to study of GT transitions to highly 
excited states. One of the advantages of the present method is that one can describe various final states in daughter 
nucleus(Be) such as low-lying cluster states and high-lying non-cluster states. Although the present model space is 
not a complete basis, it exhausts the GT transition strengths exactly. 

The present framework is a kind of bound state approximations and continuum states are not incorporated. The 
coupling with continuum states should be taken into account carefully in discussion of widths of excited states above 
the thresholds. 
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